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Abstract 

Our goal is to clarify the relation between entanglement and cor- 
relation energy in a bipartite system with infinite dimensional Hilbert 
space. To this aim we consider the completely solvable Moshinsky's 
model of two linearly coupled harmonic oscillators. Also for small val- 
ues of the couplings the entanglement of the ground state is nonlinearly 
related to the correlation energy, involving logarithmic or algebraic 
corrections. Then, looking for witness observables of the entangle- 
ment, we show how to give a physical interpretation of the correla- 
tion energy. In particular, we have proven that there exists a set of 
separable states, continuously connected with the Hartree-Fock state, 
which may have a larger overlap with the exact ground state, but also 
a larger energy expectation value. In this sense, the correlation en- 
ergy provides an entanglement gap, i.e. an energy scale, under which 
measurements performed on the 1-particle harmonic sub-system can 
discriminate the ground state from any other separated state of the 
system. However, in order to verify the generality of the procedure, 
we have compared the energy distribution cumulants for the 1-particle 
harmonic sub-system of the Moshinsky's model with the case of a cou- 
pling with a damping Ohmic bath at temperature. 
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1 Introduction 



The concept of entanglement has been recently considered by many authors 
in connection with several properties of the quantum systems and as a po- 
tential resource in quantum computation and information processing both in 
discrete and in continuous variable systems PQ[2]. Moreover, entanglement 
has also been recognized to play an important role in the study of many 
particle systems [3] and experimental and theoretical studies have demon- 
strated that it can affect thermodynamical properties both of the quantum 
phase transitions in the condensed matter and in molecular systems jl] [5] 
[BJ. However, most of the studies on the subject consider only finite di- 
mensional Hilbert models, which is not the typical situation occurring in 
atomic/molecular physics. As pointed out in [2J, the theory of the entan- 
glement for the infinite dimensional setting is full of difficulties, which can 
be cured choosing suitable subsets of the density matrices. In particular, 
the von Neumann entropy is not a continuous function in the Hilbert space, 
and for any given state of finite entanglement, one can find at least another 
state closer as we want to the previous one in trace norm, which is infinitely 
entangled. 

However, a new area of research has been opened by [6j [7j [8], where it 
was shown that the entanglement, even if it is not a quantum observable, 
can be used in evaluating the so-called correlation energy: that is the dif- 
ference between the true eigenvalue energy of a given composite system of 
identical entities, with respect to that one prescribed by the Hartree - Fock 
(HF) method. In [9] the case of the formation of the Hydrogen molecule was 
discussed, and a qualitative agreement between the von Neumann entropy 
of either atom ( as a measure of the entanglement of formation of the whole 
system) and the correlation energy as functions of the inter-nuclear distance 
was shown. However, the extension of this idea to multi-atomic molecules 
and its effectiveness remains still unclear [8] [9] . Actually, the correlation en- 
ergy is an artifact of the approximation procedure, then it is not a physical 
observable and, by second, it can be modified by the adopted method of 
calculations. Nevertheless, since any other approximating disentangled state 
has also a larger energy expectation value with respect to the HF state, one is 
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lead to look at the correlation energy as an entanglement gap in the sense of 
[7] . Adding any further correction term at the wave function has to decrease 
the energy expectation value to the Hamiltonian eigenvalue and increase the 
entanglement at the same time up to fill the gap, describing in such a way 
a peculiar domain in state space, around the exact one. Thus, the first aim 
is to quantify such a kind of relation, at least on a specific model, finding a 
quantitative expression of the entanglement in terms of the correlation en- 
ergy of the ground state for a composite bipartite system. In order to have an 
analytically tractable example containing all the desired features, we treated 
with a family of two coupled harmonic oscillators [10] in 3 dimensions. The 
coupling constant of the two parts is in a one-to-one correspondence with the 
correlation energy and and with entanglement estimators. Viceversa, assum- 
ing that such a relation is invertible, and then an estimation of the correlation 
energy can be expressed in terms of the entanglement, one may ask how direct 
measurements on one of component subsystems can provide such quantities. 
To this aim we have found for the considered model an expression of the 
concurrence, in terms of the momenta of the 1-particle subsystem energy 
probability distribution, from the knowledge that the composite system is 
in its ground state. Thus we have an entanglement witness and an a priori 
estimation of the correlation energy. However, in concrete we have to be 
able to distinguish the energy probability distribution of the entangled state 
from any other yielded by a mixed state or a pure thermal state. First we 
give an upper bound to the environmental temperature, over which all our 
procedure losses validity. Then, we compared the distribution generated by 
coupling a single harmonic oscillator with an Ohmic bath at O-Temperature, 
by the analysis of all energy distribution cumulants. 

In Sec. 2 we briefly review the main properties of the model: its exact 
fundamental state, the HF approximation and the correlation energy. Sec. 3 
for the fundamental state of the Moshinsky's model we evaluate the von Neu- 
mann entropy and the concurrence in terms of the correlation energy. Since 
the concurrence can be expressed in terms of the dispersions of observable 
conjugated quantities, also the correlation energy takes a well defined phys- 
ical meaning. A similar relation can be established for the fidelity. In Sec. 4 
we prove the existence of a continuous manifold of pure separable states, con- 
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taining the HF state, having an overlap with the exact ground state, which 
can be larger with respect to the former. In Section 5 we describe how the 
1-particle energy probability distribution for the exact state of the model can 
be distinguished from that one for a single harmonic oscillator coupled with 
an Ohmic bath at O-Temperature. Some final remarks are addressed in the 
Conclusions. 

2 The Moshinsky's Model 

In order to evaluate how good the HF mean field method is in computing 
quantum states, Moshinsky proposed a simple, but non trivial, model of two 
coupled spin-| harmonic oscillators in 3 dimensions in [TU]. In dimensionless 
unities, the Hamiltonian of the model reads 

H = \ (ft + ft + f? + ft) - \K (rl - f 2 ) ' , (1) 

where r*j and f>i denote the position and the momentum operators of the i-th 
particle, respectively. The constant K parametrizes the interaction strength 
of a supplementary quadratic potential between the two oscillators (notice 
the difference of sign with respect to [10J ) . The model describes a system 
of two identical particles in the same harmonic potential, interacting by a 
smooth effective repulsive coupling, which is truncated at the second order 
in a Taylor expansion for small interparticle distances. Thus, we will dwell 
upon < K < 1/2, where the upper bound will correspond to a breaking of 
the model, since no bound states can exist. This signals that the model is 
far to be realistic and it is intended only as a toy model shaped to our aims. 
The model energy spectrum is 

E n , m = ~ (1 + X 2 ) +n + m X 2 , m,«6NU{0} (2) 

where 

* = (l-2tf)*, (3) 

plays the role of effective coupling constant, parametrizing a sort of double 
well potential, with an increasingly higher (or wider) barrier for \ — > 0. The 
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normalized position wave function of the fundamental level is given by 



= {%y'~e-*^-~** /2 > ( 4 ) 
where the mean and relative positions 

have been defined, respectively. In general the spectrum shows degeneracies, 
but the lowest level is always simple, except for x — 0, i.e. for the limiting 
value of the coupling K — > 1/2. Moreover, crossings occur for higher eigen- 
values at isolated points of K, but we are not interested to them. Finally, 
since the function (j4]) is symmetric in the interchange r\ -<->• f; 2 , the total spin 
must be necessarily into the singlet state. Thus, the spinorial aspect of the 
problem is not relevant at this stage, and it can be ignored. 

Applying the standard HF mean field approximation for the ground state 
of the Hamiltonian ([Q), one is led to the wave function 



2„2 , 



*) = tt- 3 / 2 (1 - Kf'e-^ 1/2 { R2 ^ 2 )l\ (6) 
corresponding to the approximated eigenvalue 

Ehf = 3(1-K) 1 / 2 . (7) 
Defining the correlation energy ( positive, by Ritz's theorem ) as 

£co rr = £hf - £o,o = 3VT^K - - (l + Vi^2K) . (8) 

Moreover, the explicit expression of the overlap (or the squared fidelity) 
between the exact and the HF wave function is 

64(1 -OT- 2*)*' <lf (9) 

((l + \f\^~K) (l + \/T~ r 2K) - K) _ 

Thus, one can figure out that adding to the HF state further corrections, 
surely the estimation of the energy eigenvalue improves and the fidelity in- 
creases, but the simplest factorized expression in ([6]) will be lost. Differently 
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to what happens in the approximated state, the two oscillators in the correct 
fundamental state are entangled. From analytic point of view this happens 
because of the different coefficients of R and r in the wave function (J3J). 
From a different point of view, one can see the expressions (J1J) and (jH]) as 
two distinct continuous curves in the Hilbert space, parametrized by K (or 
X). They have only one common point at K = 0. The main property of the 
latter curve is to contain only factorized states. 

3 Entanglement Estimation 

Since we are dealing with pure states, the entanglement estimator is the 
entanglement entropy, given in term of the von Neumann entropy 

SVn [pr] = -Tr [p r log 2 p x \ = - ^ A*i lo S2 (Mi) > ( 10 ) 

i 

of the reduced to 1-particle density matrix p r = Tr 2 [p] = Tr x [p], denoting 
by fii the corresponding eigenvalues. 

On the other hand, the von Neumann Entropy S v n satisfies the additive 
relation SVn IP ® = ^vN [p] + S v n [&} , for any factorized density operator 
p®a. But this is precisely the structure of the reduced density matrix, which 
factorizes into positional and spinorial contribution, where the latter takes 
the form a — 1 1 for the singlet spin state. Thus, it contributes to an addi- 
tive constant term (equal to 1), which measures only the equal uncertainty 
in attributing one of the two possible quantum states to each spin. Following 
the ideas in [11] for fermions, anti-symmetrizing the product of 1-particle 
orthogonal states into a spin stationary state contains all information about 
entanglement by definition. In conclusion, here we will compare only the con- 
tributions to the entanglement coming from the space configurations factor 
of the 2-particle fundamental state. 

In the position representation the exact 2-particles density matrix p for 
the fundamental state (jl]) is given by the integral kernel 

Po^rVI) = (|)V(* 2 +*' 2 )/V*V^)/ 2) 
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where the supplementary variables R' and r* are in analogy with (J5J). A simi- 
lar expression holds for the |^hf) state fl6]), where the density matrix is given 
by the product of gaussian normal distributions with the same variance. The 
consequences of such different structure can be seen also by the comparison 
of the 1-particle space distribution densities, which are given by 



Poi (r) 



2X 2 



(X 2 + 1) 



7T 



3/2 



-2- 



Phfi (r) 



;i - Kf^ 



71 



3/2 



(12) 



Thus, because of the repulsive interaction, the exact average distance between 
the particles is larger than in the approximated estimation, being their ratio 



2V1-2A" 



1/2 



, with a divergence for K — > \ (x — > 0). 



The exact 1-particle integral operator density matrix p r = J p r 



r, r 



dr 1 has the kernel 



Pr {r, r) 



2X 2 



7r( X 2 + l) 



3/2 



exp 



2 (x 2 - 1) r ■ r 1 - ( 4 X 2 + (x 2 + 1) ) (r 2 + r' 2 ) 



(13) 



That can be rewritten in the usual gaussian form [1] [T5] 

*2 



p (A P ,A,) (f y 



1 



27rAg 



3/2 

, exp 




Ap {f— r 



, (14) 



where the squared mean values for the 1-particle position and momentum 

1 



X 2 + l 

4 X 2 



Ap 



4 



(x 2 + i) 



(15) 



have been introduced. 

The system becomes disentangled when the 1-particle state is the pure 
minimal packet, i.e. when Ap 2 Aq 2 = \. But this occurs only for K = (x = 
1). In this sense Ap 2 and Aq 2 contains information about the entanglement 
of the system, as remarked in jl]. However, differently from jl], by the simple 
algebraic relation x 2 = ^2 the squared mean values contain also information 
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about the form and the strength of the interaction. Thus, one would arise 
the question if the analysis of ( I14p not only provides information about an 
entangled harmonic oscillator, but also the main properties of the coupling: 
is it coupled to a small system or to a thermodynamic bath? 

From expression ({TBI the kernel of the eigenvalue equation for p r is sym- 
metric and of Hilbert-Schmidt type, since the coefficient of r 2 (and r ' 2 ) is 
negative. So the spectrum is real and discrete, as for the tensor product of 
three independent oscillators. Accordingly, the eigenf unctions of p r can be 
factorized in the product of three functions, each of them depending only on 
one real variable and the eigenvalues product 

Ul,m,n (r) = Wi{x)w m {y)w n {z), pi^ n = ViV m V w (16) 

so that the problem is reduced to solve the 1-dimensional integral spectral 
problem 

J p r (x,x') wi(x') dx = UiWi(x), (17) 
which has the non degenerate spectrum and eigen-solutions of the form 

vx = Cc l , Wl (x) = H { {^xx) e'i x \ (18) 

where C = ; c = (iTx) anc ^ {Vx x ) denote the Hermite polyno- 

mials. Of course, these positive eigenvalues sum up to 1, because they are 
related to a matrix density operator. On the other hand, accordingly to ( TT6l) 
the eigenvalues of the one-particle density matrix p r are given by 

^ = C 3 c l+m+n = C 3 c k^ fc £ ^ Qgj 

with degeneration order deg [pk] = k (k + 3) /2 + 1 . Thus, the 1-particle 
density p % (r) = J2i mn C 3 c l+m+n \wi (x) w m (y) w n (z) | 2 is represented as a 
mixed state in the basis of the "natural orbitals", using the terminology by 
Lodwin [12J, which describe a 3D single harmonic oscillator of frequency \. 
The weights pu describe a system at the equilibrium temperature k B T* = 
3x/(2hiYr^), which is a decreasing function of x- Since experiments are al- 
ways performed at a finite temperature, performing particle position/momentum 
measurements on such a system we need to work at T < T* (x), in order to 
highlight the quantum behavior discussed here. 
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For comparison, the spectrum of the reduced density matrix in the the 
HF approximation is given by {1 (simple), (infinitely many degenerate)}, 
with the corresponding eigenf unctions (for each of the three space variables) 
{Hi ((1 - Kf'x) exp [-ii^a; 2 ] }. 

Hence, if we are allowed to look at the eigenvalues of the density matrix 
operator p r as the probabilities to find the 1-particle subsystem in one of 
the states of a i^-parametrized family of harmonic oscillators, for small K 
it can be found very likely in the fundamental one. But this probability 
decreases rapidly to for K — > |, while the higher states become significantly 
more accessible. Notice that at K = A the system is meaningless, since all 
eigenvalues of p r become except one. However, for < K < | one can 
analytically sum up Tr (p r ) = 1 pointlike, taking into account the degeneracy. 
On the other hand, the lack of coherence can be estimated also by computing 
the Tr[pl\, which is 1 only for pure states. In the present case one has 
explicitly 

^ [Pi] = -f^s, (20) 
(1 + XV 

which is a monotonically decreasing from 1 to function on K. 

Complementary to this quantity there is the so-called linear entropy [T3] . 
analogous to the concurrence [1] 



C[p r .] = l-Tr[p r 2 ], (21) 

which takes values in the range [0,1]. It is invariant under local unitary 
transformations on the separate oscillators (reduced to changes of phases). 
This entanglement estimator seems to be quite useful in the present case, 
since it is always bounded, even if it is defined on an oo-dimensional Hilbert 
space. 

The entropy of entanglement ( flOl) can be explicitly written as 

SvN [p r ] = [(x + I) 2 In (x + 1) - 2x hi4 X - ( X - l) 2 hi \ X - 1|] . (22) 

For K — > | (x — > 0), the entropy S v ^ [p r ] is logarithmic divergent, according 
to the expansion 

5 vN [^]«-^jM + 0(l). (23) 
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This is a well known result for harmonic chains [T5] , indicating the degeneracy 
of the ground state in the considered limit. 

On the other hand, the behavior of 5Vn [pr] near K — >■ can be described 
by its series expansion 

SvN [Pr] = ~ 3K l { l^ K) HK) + 0(K*), (24) 

approaching 0, because of the oscillators decoupling. However, this approx- 
imation becomes inaccurate very rapidly. From the above expression, for 
K — > 0, the asymptotic behavior of the entropy is controlled by a logarith- 
mic term, differently from the correlation energy (JBJ), which has a pure power 
expansion. Then, we cannot expect a great similarity between the two func- 
tions, also at very small values of K. This results breaks the conjectured 
existence of a simple relationship between the two quantities. On the other 
hand, let us observe that both functions SVnIA] an d E COTT are monotonically 
increasing in K. Thus, the entanglement is an increasing function of the 
correlation energy (see Fig. [1]). In order to have analytic expressions, we 
solve algebraically the coupling constant in terms of E corT as 



X (-Ecorr) 



2-E'corr (yE COTT + 3) (2,E COTT -\- 3) 3E COTT (£< corr -|- 3) 



(25) 

and replacing into SVnI/v], we obtain a one-to-one correspondence S v t<s(E cott ). 

In particular, one can look for asymptotic expressions of the entanglement 
for small values of -E cor r; corresponding to small couplings. Indeed, including 
logarithmic corrections at the lowest order near E corr w 0, one obtains 

5vN (-Ecorr) ~ ^ln(2) " ^ corr + ( E clri) > ( 26 ) 

for the Moshinsky's oscillators. One verifies that similar expressions can 
be obtained studying other systems (for instance the 2-points Ising model), 
but up to now does not exist a general procedure to compute directly the 
coefficients appearing in the above developments. Moving to the upper limit 
K —t 1/2 (x — > 0), or equivalently -E C orr - > E covv = | (— 1 + \/2), the entropy 
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Figure 1: The entanglement as a function of the correlation energy for the 
Moshinsky's model. 



diverges logarithmic as 

S vN (E COTT ) « - 31n(jSc 1 ° ri rr ( ~ EcOVv) +0(1), (27) 
which is the specific behavior for degenerate ground state, as remarked for 

(E3D- 

However, the singular behavior near of the entanglement as a function 
of the correlation energy does not seem related to the specific way of its esti- 
mation. In fact, by using Eq. ff25|) into (120]) . as a function of the correlation 
energy the concurrence for the Moshinsky's model takes the form 

c < £ -> = 1 -p^ < 28 > 

/ = 9 + 12 E co „ (E co „ + 3) —4 y/2 J E co „ (E COTT + 3) (2_E corr + 3) 2 , 

a graph of which is shown in Fig. [2j This function is regular in the origin, 
but it is not in its second derivative. Again a singularity is signaling a 
faster increase of the entanglement for small values of the correlation energy. 
But, the expression fl28l) is algebraic and it can be manipulated more easily. 
Specifically, for small values of the concurrence one gets the correlation energy 
as an half-integer power series of C 

E cott + y|c 3 / 2 + ^ + 0(C 5 / 2 ), (29) 
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Figure 2: The entanglement expressed in terms of the concurrence as a func- 
tion of the correlation energy for the Moshinsky's model. 



while for <C C < 1 the expansion is 



E corr « E COTT - - C) 2 / 3 + 1 (V2 - 4) (1 - C) 4 / 3 + O ((1 - C)) 2 . (30) 



These expressions give direct relations between the correlation energy 
and the entanglement, which may suggest an experimental measure of the 
entanglement and of the correlation energy. In fact, let us suppose to perform 
two independent series of measurements of position and linear momentum on 
one particle of the system. Their results are distributed with squared mean 
values Aq 2 and Ap 2 , respectively. On the other hand, by resorting to the 
relations f JT5|) in terms of the coupling constant \ and to the expression (120]) , 
one gets 



Thus the entanglement is related to the ratio between the uncertainty and 
the energy mean value of the observed subsystem. Furthermore, by the 
definition (1211) and in the range of validity of expansion (129]) (or (1301) ), one 
may obtain a relation among observable quantities and the mathematical 
artifact -E CO rr- On the other hand, relation ( 1311) has to be used carefully 
since, if applied to a generic gaussian separate pure state, it does not give 
a measure of entanglement, of course. The point to be remarked is that its 




(31) 
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value depends by the special relation of the mean squared deviations on the 
coupling constant, not only on the preparation of the state. 

Finally, the fidelity of the fundamental state of the Moshinsky's model 
with the corresponding HF state or, equivalently, the overlap can be 
expressed as a function of the entanglement. In some sense, we are comparing 
two different ways to measure the "distance" between the two curves of states, 
even if neither quantities actually have the properties of a distance. However, 
also in this case a monotonic function can be obtained for any pair of states 
corresponding to the same coupling constant K, or correlation energy -E corr . 
This function is regular, even if at the extrema a singularity in its higher 
derivatives appears. 

4 The gap of entanglement 

Here we would like to elucidate the special role played by the HF state 
in the set of all separable states, which may be closer, in the sense of the 
trace - norm, to the exact solution. To this aim and since we are looking 
to a neighborhood of the ground state in the Hilbert space, let us restrict 
ourselves to the pure separable states, which are symmetric with respect to 
the change lf>2 (the spins are into the singlet configuration) and generated 
by wavefunctions of the form 



where for convenience we assume that $ is normalized to 1. Of course, more 
general choices are possible, compatibly with the assumed identity of the 
particles. In the class of states ( 1321) there exists the 1-parameter curve given 
by the gaussian functions 



certainly containing the HF wave function fl6]). Its overlap with the ground 
state (HJ is such that 




(32) 




(33) 



64a 3 x 3 



> \(*hf\*o)\ 2 



(34) 



(a + l) 3 (x 2 + a) 



3 
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(see eq. ©) for a G 



Of course, the 



a^(x),y / |(x 4 -l) + l = VT^K 

upper bound is exactly the value involved in the HF wavefuction (see eq. 
fl6])). The lower bound cti ow (x) is an algebraic positive monotonic increasing 
function, going from to % — > 0, like a>i ow (x) w 7% 2 + [x 4 ] ( 7 = const) 
, to 1 for % 1 (i^ 0), when the two extrema coincide. In particular, 
one sees that the maximum of overlapping is achieved for a max = X-> f° r 
which one has \(^ am J^o)\ 2 = pg^e, equal to 1 only at X = 1 (K = 0). 
It should be noticed that a max is exactly the same exponent appearing in 
the eigenfunctions of the 1-particle reduced density matrix operator (see eq. 
(j!8p ). accordingly with the notion of "natural orbital" . In conclusion, the HF 
state is not the closest (in the sense of the trace-norm) pure separable state 
to the exact fundamental state and one may wonder if other states arbitrarily 
close to it may be found. Of course, by deforming the pure gaussian form 
fl33l) with maximal overlapping, in the base of the Hermite polynomials one 
can construct symmetric factorized wavefunctions of the form 




*c = n N(i) ( E c T H i {vx4 ] ) ) ( E 4 )H j (v^?) ) (35) 

o=o 



for arbitrary complex constants | c j^| an d for suitable normalization con- 
stants iVW. Thus, it is not difficult to find 



The three factors in the r.h.s. of the above expression are < 1, then the 
overlapping of the generalized wave-packets (jHSJ) with the exact state cannot 
exceed the maximal one. In conclusion, we have proved that there exists 
a dense set of pure and separated states, containing the HF state, having 
overlap \(^ HF \^ )\ 2 < |(^|^ )| 2 < 1, except for x = 1 (K = 0), when 
^ hf = ^o- Then, the exact state cannot be approached arbitrarily close by 
a separated state, except when it is itself separate. This result is complemen- 
tary to the statement that entanglement entropy of a continuous model is a a 
discontinuous function, diverging at infinity in any neighborhood of any pure 
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state [T5] • The maximal overlapping is provided by taking the suitable tensor 
product of the natural orbitals [12]. Finally, because of the convexity of the 
set of all separable mixed states, i. e. of the form p = J2 n Pn\^c„)(^c„\, with 
p n > and J2 n Pn = 1) one can extend the previous statement to the entire 
space of states. 

On the other hand, the HF state has been selected in the class of sepa- 
rable states by the minimal energy requirement. But in the domain of the 

pure separated states of form (1331) the relation (H)^, a — E HF = — — ^- — — 
holds. Then, the expectation value of the Moshinsky's Hamiltonian H gets 
its absolute minimum indeed at the HF state. Moreover, this can be seen 
also considering general factorized states as in f l35l) . Now, because of the 
convexity of the set of separable states, the minimum in the spectrum of a 
bounded observable from below is always achieved by a pure separable state. 
Thus, we conclude that the above introduced correlation energy is not simply 
a mathematical artifact, but it looks analogous to the concept of entangle- 
ment gap introduced in Ref. [7]. Since this is a global result, not depending 
on a particular computation procedure, we claim that the HF state for the 
Moshinsky's model provides the minimum separable energy as introduced 



in 



[7] E sep = min ahle Tr Hp 



'separable 



Ehf- Moreover, the observable 



Zew = H — Ehf^- is an entanglement witness, the spectrum of which is non 
negative on all separable states and there exists the ground state (entangled) 
of the Moshinsky's model for which its expectation value is —E corr < 0. Ac- 
tually, for x > (K < 1/2) isolated eigenvalues of Zew m &y exist in the gap 
[—E corr ,0[, but they corresponds to higher energy entangled states. Thus if 
in a measurement of Zew we obtain a negative value, we can still affirm that 
the the system is in an entangled state, even if not necessarily in the ground 
state. 



5 Entanglement Energetics 

In the previous Sections we have shown that the entanglement gap is the 
main energetic scale that dictate if a composed system is, or not, entangled. 
For the Moshinsky's model we have shown that this gap is given by the 
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correlation energy, derived from the HF calculations. However, this relies 
on the full knowledge of the density matrix, while for pure states all needed 
information is encoded into the reduced matrix of a selected subsystem: in 
our case one of the harmonic oscillators. Thus the question if one can estimate 
the entanglement by energy measurements on the single harmonic oscillators 
arises, conditionally to the knowledge that the whole system is not in a 
separated ground state. These will be subjected to statistical fluctuations, 
which in principle contains the required information, i.e. the entanglement 
of original ground state of the composite system. This Section is devoted to 
how extract this result and how to distinguish the energy distribution of the 
entangled composite system, from the effects of couplings to a more generic 
environments, like an Ohmic bath, even if the latter is at temperature. 

First step concerns the calculation of the energy distribution for one single 
harmonic oscillator included into the Moshinsky's model. To this aim it 
is useful to have the above expressions in the simple harmonic oscillator 
Hamiltonian eigenvector basis B = {v/ mn } ( i — 1,2). 

First, the overlap of the exact wave function with a generic factorized 
state can be evaluated from the set of amplitudes 

(...1/2 \ 3 J T J 
A \ -'£,/ / -'m,m / -'n,re / (Q7\ 
vr J V2^ 2 ^2^'/!/'!m!m'!n!n'!' ( } 

where the matrix {I m ,m'} has a sort of chessboard structure, given by the 
relation 

I m , m> = 2ire (to + to') {-if (to + to' - 1)!! (1 + (f 2 C^, (38) 

where the expression ( = and the scaled step function e (to) = | ^ odcT 

have been introduced. 

On the other hand, in the basis B the amplitudes of \l/ a are given by 

/I 2 I T \ Oj ^ A[ Ajj-^Ayyjf A n A n f , , 

(VW^Vn>'l a) = ^-J========= : (39) 



where 



A l = e (l)2 l+1 ^a- 1 - 1 (2 - a 2 f 2 T ^-±^j (40) 
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with a = a/1 + a for brevity. The expression f)40p establishes that the non 
vanishing terms occur only for even principal quantum numbers 1,1', ... , but 
they are not correlated among them. 

In the representation (|37p - (|38p the elements of the full matrix density 
operator for the exact ground state are given by taking the tensor product 
in three dimensions of the 1-dimensional factors 

l + l' + m + m' 

Pl,V,m,m> = 4(1+C)XC 5 (41) 



d'+m' 



e(l + V) e(m + m') (/ + /'- 1)!! (m + m' - 1)!! 



[2 l + l '+ m + m 'l\l'\m\m'\] 1/2 

The corresponding reduced density matrix p r ^ v = ^2 m Pi, m ,i',m can be com- 
puted from the above expression, or using the continuous basis representa- 
tion, contracting with respect the suitable states of the uncoupled harmonic 
oscillator. In particular, we are interested in the evaluation of the diagonal 
elements p r ,i,i, which represent the probabilities to find the 1-particle sub- 
system into the energy eigenstates of the uncoupled harmonic oscillator. It 
results that these quantities are related by the following recursion relation 

^ + M +1 =(^x(C + l)C'^(^ y?7 ) (42) 
with the expressions for the fundamental state 

"*° = 2 ^f=| ■ (43) 

Thus, the general structure of the considered distribution is 

Pr,i,l = 2C l X iC+ ll Qi (C) , (44) 
(4-C 2 )^ 

where Qi (() is a polynomial of degree I in the "scaled" coupling constant 
(. This distribution of probability has its own peculiarities, which make it 
different from a generic factorized state or from a pure equilibrium ther- 
modynamical distribution. Then, for comparison one computes the energy 
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probability distribution for a factorized gaussian state by using the formulae 
(|39|) - fj40l) . For each eigen-state label one obtains the expression 



1/2 e (l) 2 (1 - 1)\\ 2 fl-a\ l , N 

p rll = 16a 1/2 ^— i ^— . 45 

The first observation is that this distribution is different from only for even 
I: this is a common character of all factorized gaussian states, included the 
HF approximated wave-packet, so it could be used to make an experimental 
comparison with the entangled state. 

On the other hand, one can compute such a kind of quantity by the direct 
use of the 1-particle reduced matrix (|T4|) [4j. In fact, by using the generating 
matrix for the Hermite polynomial, the diagonal elements of the reduced 
matrix in the basis of the pure harmonic oscillator are given by 



PrU " [(2Ap 2 + 1) (2Ag 2 + l)f +1 ^ 1 \ v /(4Ap 4 - 1) (4Ag 4 - 1) 



(A P ,A 3 ) _ o t( 2A P 2 - 1) ( 2A ? 2 - l )f 2 P ( 4Ap 2 Ag 2 
' " ■ ' -'-'--"''I 

(46) 

where Pi denotes the Z-th Legendre polynomial. The parameters Ap and Ag 
are independent quantities, limited only by the minimal uncertainty condi- 
tion Ap 2 Aq 2 > \. Of course, substituting the expressions of Ap and of Aq 
given in ( fT5|) , one recovers the formula ( 144]) : there the emphasis is on the de- 
pendency by the coupling strength. In Fig. ([3]) we give a set of contour plots 
of the probabilities to find the single harmonic subsystem ( of frequency 
u) = 1) in one of the first six eigenvalues as functions of the uncertainties 
(Ap, Ag), accordingly to expression ( ]4"6]) . The bold dashed curve is given by 
the equations (|15|) of the uncertainties in the Moshinsky's model. Of course, 
the efficacy of the above procedure to measure the entanglement has to be 
evaluated by comparison with other situations. For instance, one may ask if 
is it possible to distinguish the above distribution of energy eigenvalues from 
a sufficiently general mixed one. Specifically we consider that one obtained 
coupling one of the harmonic oscillator to an Ohmic bath [16]. To this aim, 
we propose two different methods for this comparison. 

The first way is based on the position and momentum measurements, from 
which we can reproduce the graph of Fig. El for fixed values of the coupling 
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Figure 3: The probability distributions fj4*6]) to find into one of the first 
six modes of the simple harmonic oscillator a 1-particle subsystem, de- 
scribed by a reduced density matrix ( TT4"j) in the (Aq, Ap) plane. In the 
area between two contours the probability varies by 1/10 of the maxima 
values {0.942809, 0.249761, 0.190042, 0.105418, 0.103428, 0.0669486}, respec- 
tively for each plot, decreasing from black to lightest gray. The black thick 
dashed curve represents the values of (Ag, Ap) given by f|T5|) . while the black 
thick continuous curve corresponds to the oscillator coupled to an Ohmic 
bath at T = 0° accordingly to (jSO^gThe dashed gray boundary curve is 
determined by the hyperbola AqAp = ~, which represents the gaussian fac- 
torized states. When the coupling constants of both models vanishes, the 
corresponding wave-packets are minimal. 



constant. In this case, the curves corresponding to the parametric represen- 
tations of (Aq, Ap) characterize the two models. Thus, we can distinguish 
the Ohmic model from the Moshinsky's model by knowing the position and 
momentum uncertainty behaviour in the (Aq, Ap) plane. Because of such 
measurements, this method produces a lack of information about the ener- 
getics of the system. Then, the second method, suggested by [I] , is based on 
the analysis of the cumulants of the simple harmonic oscillator Hamiltonian 
H HO , namely 

m n HO )) = (-ir—\nz(i) 

where the partition function Z(£) is evaluated by tracing the harmonic prop- 
agator in the imaginary time £, with respect the generic gaussian density 
matrix ( fl4l) . Following the standard calculations [I] [17], one knows that 



(47) 



<H HO \ 

Pr 

3/2 

(Ap 2 + Aq 2 ) sinh£ + 2Ag 2 Ap 2 (cosh£ - 1) 



Z(0 = (e-^°) o (A P ,A,) = (48) 

1 + cosh £ 
2 _ 

Then, a list of the cumulants can be algorithmically computed as polynomials 
of even dergee on the uncertainties (Aq, Ap), as for instance 

((H% )) = \{(n- 1)! (Aq 2n + Ap 2n ) + O (Aq 2 ^ + Ap 2 ^)) , (49) 



For the Moshinsky's model one substitutes the uncertainties given in (1151) 
in (j49p . while for the Ohmic bath (in the underdamped limit) one uses the 
mean squared values jl] 

v ,- n 2 ' 



Ap 2 = (1 - 2(3 2 ) Aq 2 + W In ^, (50) 

7T UJ 

where ft is the coupling to the dissipative environment, in units of the oscil- 
lator frequency, and uc is a cutoff frequency. 

In order to have a unique parameter, which measures the interaction 
strength between the singled out oscillator with the remaining of the com- 
posite system (the environment) in both considered cases, let us assume the 
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Figure 4: The second energy cumulants {{Hfj )) as a function of the coupling 
constant \ with uc = 10 u. Dashed line represents the Ohmic bath, while 
solid line corresponds to system ([TJ. 



This relation is suggested by the derivation of the classical Langevin equation 
from a pure Hamiltonian system of coupled oscillators [16J. Thus, we can 
rewrite the above cumulants in terms of the coupling constant \ for both 
systems and compare them, as it has been shown in Fig. H] for those of 
order 2. There, we can see that the two functions are similar for x = 1, 
i.e. when they describe free harmonic oscillator in both cases. The situation 
dramatically changes when \ decreases, i.e. for stronger interactions. In 
fact, for x ~~ >* 0, the second energy cumulant associated with the Ohmic 
bath remains finite, while for the Moshinsky's model it is divergent. Similar 
considerations can be made for the higher order cumulants. Thus, we have 
provided a method for distinguishing the two classes of states. 

A different approach concerns the analysis the logarithm of the cumulants, 
at a fixed value of the coupling parameter x, as function of their order n. In 
fact, it results that \n((Hjj )) is approximatively a linear function of n. But, 
from expression ( 149|) . the relevant physical information is contained in their 
slope and in the corresponding differences between the two models. Then, 



relation 



(i-x 4 ) 2 

4(1 + X 4 )' 



(51) 
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Figure 5: The relative slope difference R (x) as a function of the coupling 
constant X- 



we introduce the function 

which gives the relative difference of the \n{(H^ _ 7 )) (J = Oh, Mosh) slopes 
for the two models at different x values (see Fig. \5§. By inspection, we can 
deduce that for stronger or weaker interactions, i.e. for x ~ or X ~ 1> 
the two models are well distinguishable. While this becomes less obvious for 
an intermediate range of the coupling constant, where the relative difference 
takes values R(x) < 0.5. 

Finally, because of the explicit dependency of the cumulants from the 
position/momentum uncertainties, one may express the latter in terms of 
the first (mean energy) and second (variance) energy momentum. Thus, 
adopting the formula ([31]) as a common expression for any model of harmonic 
oscillator coupled to an environment, one provides a new expression of Tr [pf\ 
as 

N " w^hJ) ) (53) 

and then of the concurrence only in terms of measured energy distribution 
properties. Moreover, from the arguments of Section 3 we get also an a priori 
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estimation of the correlation energy. Now, using the parametrization of the 
coupling constant (|51|) . one can compare the resulting concurrencies for the 
considered models (see Fig. |6]). Qualitatively one can establish the intensity 
of the entanglement generated in those different models. 

6 Conclusions 

In the present article we have clarify the relation between the entanglement 
and correlation energy in a bipartite system with infinite dimensional Hilbert 
space. We have considered the completely solvable Moshinsky's model of two 
linearly coupled harmonic oscillators, which may constitute a simple case be- 
fore studying more complicated systems, like double well potentials. The 
system has a coupling constant, which can be varied in a finite range. Thus, 
it continuously parametrizes two special curves in the space states: one con- 
taining the exact ground states, the other the separable HF states. Of course, 
for vanishing coupling the two curves emerge from the same state, but their 
separation can be described in terms of norm, entanglement and energy cor- 
relation. The peculiarity of the second curve is to lie always in a set of 
entropy entanglement states, while along the first one it increases monoton- 
ically (in K), with a logarithmic divergence when the ground state becomes 
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degenerated. On the other hand, a similar description is given in terms of the 
correlation energy, which in principle is defined only for pairs of correspond- 
ing states (at the same coupling constant) in the two curves. We have proved 
that entanglement and correlation energy are one-to-one along these curves, 
at least for the considered model. However, they are not simply proportional, 
but at small couplings they have a quite different rate of increasing. This 
phenomenon occurs not only if one uses the entropy of entanglement for pure 
states, but also if one introduces the concurrence. However, in the considered 
model certain algebraic approximated expressions of the correlation energy 
in terms of the concurrence are given, so that an artifact of the calculation 
methods can take a physical interpretation. However, at the moment we 
have not a general method to compute directly the coefficients of such type 
of expansions. These could be very useful in order to have an alternative a 
priori estimation of the errors made in numerical computations of the correct 
expectation values of the energy. Such a type of relation may be useful in 
the studies of bipartite systems with many inner degrees of freedom, like the 
dimers of complex molecules (see [9] for instance). In this respect the ex- 
plored concept of entanglement gap and its identification we made with the 
correlation energy may play an important role: it represents the energy range 
we have to be able to measure, in order to establish if a composed system is, 
or not, entangled. However, the theory of the entanglement gap for systems 
with infinite dimensional Hilbert space does not seem completely developed 
as for the finite dimensional case and further investigations are needed. In 
the final section we have shown that, conditionally to the knowledge that the 
whole system is not in a separated ground state, one can estimate the en- 
tanglement by energy measurements on the single harmonic oscillators. This 
can be done by two sequence of position and momentum measurements, as 
well by energy measurements. The distribution of the energy measurements 
is sufficiently characterized in terms of its cumulants. This analysis enable us 
to compare among different systems at temperature and distinguish their 
ability to generate entanglement, for instance by using the parameter R in- 
troduced in ( |52i) . Finally, via the formula ( 1531) we propose a new estimator 
of the entanglement, based on the first two momenta of energy distribution 
of the considered subsystem. We plan to check the how good is the present 
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approach in considering bipartite multi-particle systems and non linearly cou- 
pled systems. In particular, we would like to consider integrable systems, like 
in [18] , in which a complete analytic control of the calculations is at the hand. 
Another direction of research is to consider a different entropy entanglement 
parameter, like the quantum version of the Tsallis entropy [T9] . 
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